Nonsingular black holes in quadratic Palatini gravity 
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We find that if general relativity is modified at the Planck scale by a Ricci-squared term, elec- 
trically charged black holes may be nonsingular. These objects concentrate their mass in a mi- 
croscopic sphere of radius r^ore ~ Ng^^lp/3, where Ip is the Planck length and A^, is the number 
of electric charges. The singularity is avoided if the mass of the object satisfies the condition 
Mq « mpQem-'V|/2, where mp is the Planck mass and Qe„i is the fine-structure constant. For as- 
trophysical black holes this amount of charge is so small that their external horizon almost coincides 
with their Schwarzschild radius. We work within a first-order (Palatini) approach. 
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INTRODUCTION 

The combination of general relativity (GR) and quan- 
tum theory leads to important results related with the 
existence of compact astrophysical objects. When nu- 
clear reactions are no longer effective in the interior of not 
too massive stars, gravitational collapse can be avoided 
thanks to the quantum degeneracy pressure of matter, 
which leads to the formation of white dwarfs and neutron 
stars pj. For very massive stars, however, this quantum 
pressure is unable to counterbalance the gravitational 
pull and matter collapses yielding a black hole. Accord- 
ing to GR, black holes are extremely simple objects that 
can be fully characterized by their mass, charge, and an- 
gular momentum and that are shielded by an event hori- 
zon, which generates the emission of quantum particles 
via Hawking radiation Hidden inside the event 

horizon, the mass and charge of the collapsed object are 
concentrated on a point (or a ring if there is angular mo- 
mentum) of zero volume. In this infinite density point 
(or ring) curvature scalars diverge, thus implying the ex- 
istence of a singularity. In such an extreme scenario, 
one would expect that the quantum properties of the 
gravitational field could halt the collapse and provide a 
new stable object in much the same way as white dwarfs 
and neutron stars arise when the quantum degeneracy 
pressure of matter dominates. On dimensional grounds, 
quantum gravitational effects should become non negli- 
gible at densities of order pp ~ c^/hG^ ss 10^^ g/cm^ [4[. 
This exceedingly large density sets the scale at which GR 
should be replaced by some improved description able to 
stabilize the collapsed object. Due to the lack of a fully 
developed and tractable quantum theory of gravity, we 
could gain some insight on how the internal structure of 
black holes is modified by considering departures from 
the dynamics of GR at high curvatures. In this sense, 
extensions of GR formulated in a first-order, or Palatini, 
approach seem a very promising framework to explore 



aspects of quantum gravity phenomenology [^, 0] • This 
type of theories have been studied in cosmological mod- 
els finding that the big bang singularity can be avoided 
in very general situations [8]. In these scenarios, the uni- 
verse began in a contracting phase that bounced off to an 
expanding phase after reaching a minimum volume char- 
acterized by a density of order p p . The successful results 
of that approach suggest that black hole interiors could 
be modified avoiding the zero- volume singularity of GR. 

In this paper we show that electrically charged black 
holes in an extension of GR defined by the Planck-scale 
corrected action {Ip = \/ fiGjt? = Planck length) 
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formulated a la Palatini, i.e. by assuming that met- 
ric and connection are independent entities [sl, concen- 
trate their mass in a compact sphere of finite volume. 
Theories of this type are expected to arise when quan- 
tum effects are considered in general curved space-times. 
For external observers, these black holes look essentially 
like those found in GR, though for a certain charge-to- 
mass ratio the geometry is nonsingular, which sheds new 
light on how black hole singularities could be removed 
by nonperturbative quantum effects. In ([1]) we have de- 
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X - r;^;,r^^, and Sm[g^.u, ^] rep- 

resents the matter action. To remain as close to GR as 
possible, the independent connection only appears in the 
gravitational sector of the theory through the definition 
of -R"0p.. 



FIELD EQUATIONS 

Let us now derive the field equations and discuss the 
peculiar dynamical properties of the theory ([T]) . Taking 
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independent variations of ([T]) with respect to metric and 
connection yields 
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v. W^ifRg^-" + 2fQR^^)] = 0, 
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where = 87r/|,/n, f = R + Pp {aR^ + Rf.^R^"'), Jr = 
Oaf, Q = Rfj^uR^'^ , and /q = dqj . Defining the tensor 
P^'^ = R^ag°"^ , one finds that ([2]) establishes an algebraic 
relation between the components of P^" and those of 
T^"^ = r^„5"", i-e., P^.'' - P^.^tJ). Once the explicit 
relation P^'^(Ta^) is known, one can express R and Q in 
terms of the matter according to the identities R — P^^ 
and g = P^"P„^ (see 1,0 for details and applications), 
which allows to solve the connection equation ([3]). The 
independent (symmetric) connection T'^^ turns out to 
be the Levi-Civita connection of an auxiliary metric h^^i, 
related with the physical metric as follows 



V = VdetE(E-i) V. 
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where Sq,"^ = fn^a + '^IqPo' is also a function of T^^ 
and E is the matrix representation of Ea''- It is im- 
portant to note that when T^,y = the field equation 
([2]) boils down exactly to GR in vacuum. This can be 
shown by taking the trace of ^ with g^^^ , which for 
our model gives R — —k^T in general and R^ac = in 
vacuum, and then using the vacuum relation P^"\vac = 

" + to obtain Q,ac = and 

Pfi^lvac — 0. From this it follows that in regions where 
Tf^v = one has l^a" = (5^, ^/^i/ = ff/^i/ and R^i, = 0. This 
shows that in our theory the dynamics departs from that 
of GR only in regions that contain source^ Moreover, 
this fact guarantees that there are no new propagating 
degrees of freedom and, therefore, the resulting theory is 
not affected by ghosts or other dynamical instabilities Q , 
which clearly manifests the inequivalence between met- 
ric and Palatini formulations for higher order gravity (see 
for instance for a discussion on this point, and [Tl| 
for a recent attempt to describe Kalb-Ramond fields in 
terms of torsion in a Palatini-Cartan approach). Thus, 
in vacuum hfj,i, — g^i, implies that F^^ coincides with the 
Levi-Civita connection of g^v However, the presence 
of matter-energy induces a relative deformation between 
/i^ii/ and g^v, given by (|3]), such that F^^ is no longer 
metric compatible. Physically this means that a given 
background metric can be deformed locally by the pres- 
ence of matter in a way that depends on how the energy- 
momentum density is distributed. This modification of 



^ In general, for Palatini f(R, Q) theories the vacuum field equa- 
tions boil down to GR with an efli^ective cosmological constant. 



the gravitational dynamics is in sharp contrast with the 
more standard metric approach, in which the connection 
is constrained a priori to be the Levi-Civita connection of 
the metric. In that case, the modified dynamics is char- 
acterized by the existence of new propagating degrees of 
freedom which manifest themselves in the form of effec- 
tive dynamical fields or higher derivatives of the metric. 
In our theory the modified gravitational dynamics is due 
entirely to the very presence of matter and its active role 
in the construction of the connection. 

The Schwarzschild black hole is the most general spher- 
ically symmetric, non-rotating vacuum solution of GR 
and also of However, that solution assumes that all 
the matter is concentrated on a point of infinite density, 
which is not consistent with the dynamics of ([2]). In fact, 
if one considers the collapsing object as described by a 
perfect fluid that behaves as radiation during the last 
stages of the collapse, explicit computation of the scalar 
Q = RfivR^" shows that the energy density p is bounded 
from above by Pmax = Pp/32, where pp = ic^ / AnhG^ 
[1, 0] • Therefore, one should study the complicated pro- 
cess of collapse of a spherical non-rotating object to 
determine how the Schwarzschild metric is modified in 
our theory. To avoid such an involved analysis, which 
would require advanced numerical methods, we consider 
instead static vacuum space-times with an electric field, 
which allows to make significant progress using analyti- 
cal methods only. Disregarding the rotational degrees of 
freedom, static spherically symmetric electrically charged 
black holes provide the simplest scenarios to address the 
qualitative and quantitative changes produced by Pala- 
tini gravities with Planck-scale corrections on the very 
internal structure of black holes. The resulting solutions 
should therefore be seen as Planck-scale modifications of 
the usual Reissner-Nordstrom solution of GR. We thus 
take 
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where F^^i/ = dfj,Ai, — d„A^ is the field strength tensor 
of the electromagnetic potential A^. We shall focus on 
purely electrostatic spherically symmetric configurations, 
for which the only non vanishing component is 
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which is a solution of the (sourceless) equations V^F'*" 
assuming the diagonal metric ds^ — gttdt^ + grrdr^ 
r^dH.^. The associated stress-energy tensor becomes 
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where / and are the 2x2 identity and zero matrices, 
respectively. To solve for the metric, it is convenient to 



rewrite © as P^.^S^'' 



^T^" and use (g]) to get 
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an equation involving h^i, and the matter only (recall 
that P^" = R^cg^n 



V det S ^ ^ 



(8) 



Inserting ([T]) into (jS]) and knowing that in our theory 



Q = 



;f , T^^" = diag[a-,a-,a+,a+], where 



a±=[l± ^ j , we find [12| 

which recovers GR when /p — 0. 

METRIC AND KRETCHSMANN SCALAR 

Since the right-hand side of (HJ does not depend on 
g^iv, we can solve directly for h^y and then use ^ to 
obtain g^^. Defining an auxiliary line element for /i^^ as 

with A{f) = 1 — 2M{f)/r, we find that ^(f) =constant 
can be absorbed into a redefinition of the time coordinate, 
like in GR. The metric g^i, can thus be expressed as (with 

z = r/^rqlp) 



gtt 



A{z) = 1 - 



A{z)_ _ a+ 

tT+ (J-A{z) 

[l + <5iG(z)] 



(10) 



X 1/2 



where we used the relation = r^(T_, and defined 



5i 




rs 



(11) 



Here = 2Afo, and Afo is an integration constant that 
coincides with the Schwarzschild mass in the vacuum 
case. From PU)) we see that all the information about 
the geometry is encapsulated in G(z), which satisfies 



dG 

dz 



(12) 



It is easy to see that for z ^ 1 we have G{z) w 
-l/z-S/lOz^ which leads to A{r) « l-rs /r + rl/2r'^ ~ 
rsrqlp/2r^ and recovers GR when r ^ Ip. Note that for 
the electromagnetic field the function dG/dz is not real 
for z < 1. This is a dynamical consequence of the theory, 
not a coordinate problem, and will be discussed later. In 
general, (|12p can be exactly solved using hypergeometric 
functions. Focusing on the small z region, which is the 
relevant one to look for departures from GR, we find 



Giz)=p + 



(13) 



where fx{z) = 2i^i[5, A, |, 1 - z% and ^ « -1.74804 is 
a constant resulting from matching the large and small z 
expansions. Expanding (|13p around z ss 1, we find 



gr 



$2 262 {S2 - (5l) 

(l + /?5i)yi^~ (l + /35i)' 



(14) 



The above expansions show that the gtt and grr compo- 
nents are in general divergent as 2: — > 1. However, they 
also point out the existence of a particular combination 
of mass and charge for which the above expansion should 
be reconsidered. In fact, if we take (5i = — 1//3 = Jjf, the 
expansion near z k, \ becomes 

gtt « (1 - ^11^2) - (1 - 25*/3J2) (z - 1) + . . .(15) 



1 1 (l~7^iV3(52) 

'2(\-6llb2){z-\) 4(l-5iV<52)' 
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Though the component in ()15p diverges as z — ^ 1, 
the introduction of a new radial coordinate z*(z) = 
2\J z — 1/(1 — ii\l&2) = 'T* I \jTqlp turns the physical line 
element near z = r j yjvqlp = 1 into 



ds^ 



(1 - 51/82) [-dt' + {dr*Y] + {rqlp)dn' . (16) 



This Minkowskian line element implies that the geometry 
at z = 1 is regular. This contrasts with the general case 
<5i 7^ 51 where, defining f = 252r/{l + fi5i), (HH) leads to 

ds^ ^ r + (rf^)2i + (^r„lp)d^f , (17) 

4(32 V 2 — 1 

which is singular aX z — 1. In fact, one finds that the 
Kretschmann scalar of blows up at z = 1 if 5i 7^ (5J, 

i?"/3p.i?a''^" - + (1 - 5i/5l)Ki + (1 - 6i/5lfK2 , 

(18) 

because the functions Ki and K2 diverge near z = 1 as 
~ l/{z — l)'^/^ and l/{z — 1)'^, respectively. Since 
i^o is finite, near z = 1 and for 5i = (5]' one finds 

R'^P^.uRJ'"' - (l6 + - ^) /r2/|,. The scalars 

g'^'^Rfii, and R^'^ R^^i, of g^i, have a similar behavior [l^. 



PHYSICAL PROPERTIES 

Let us now discuss the physical interpretation of these 
results. At the surface z = 1 or, equivalently, r^ore — 
y/vqlp, we have found a singularity that can be avoided 
if and only if the condition = 5^ is satisfied. Ex- 
pressing the charge as q = NgC, where e is the elec- 
tron charge and A^^ the number of charges, we can write 
Tq — y/2a^Nqlp, whcrc ttern is the fine structure con- 
stant. With this we find that the area of this surface 
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is given by Acore = N'qy/2a^Ap, where Ap = Anl'j, is 
Planck's area. This suggests that each charge sourcing 
the electric field has associated an elementary quantum 
of area of magnitude ^/2a^Ap. From this it follows 
that the ratio of the total charge q by the area of this 
surface gives a universal constant, pq = q/{4:Trr'^g^^) — 
{\.T^\fi)~^ ^JWJi^GP^, which up to a factor \f2 coincides 
with the Planck surface charge density. Furthermore, the 
regularity condition 8\ = 8\ sets the following mass-to- 
charge relation 

1 \r\ Mo 1 mp 

which indicates that the matter density inside the sphere 
of radius Vcore is another universal constant, p*ore — 
Mo/Vcore = Pp/4i5i, independent of q and Mq. 

These results picture a black hole interior fully com- 
patible with the discussion on compact objects presented 
in the introduction. Planck-scale physics acts so as to 
stabilize the collapsed object yielding a compact core of 
radius rcore, which contains all the mass in its interior 
and all the charge on its surface. The latter point is sup- 
ported by the linear growth of Acore with Nq and the 
fact that the auxiliary geometry h^i, is not defined in the 
region r < Vcore {z < 1). In fact, since the physical ge- 
ometry, defined by 5^1,, is completely regular at z = 1, 
the impossibility of extending /i^^ into z < 1 indicates 
a limitation of the electric field to penetrate inside that 
region. For a description of the geometry in 2; < 1 one 
should specify the T^"^ of the sources that carry the mass 
of the core and the charge that generates the external 
electric field, which would allow to define a new auxiliary 
metric h^i, able to parameterize the internal geometry 
of the core (assuming that suitable matching conditions 
can be found at z = 1). Since the regularity condition 
(IT^ and cosmological models Q support that the inte- 
rior matter density is bounded, we expect the existence 
of completely regular interior solutions. 

From the large z expansion, we saw that the GR solu- 
tion is a very good approximation for any r ^ Ip. This 
implies that the location of the external horizon of these 
charged black holes is essentially the same as in GR, i.e. 

r+ = rs (1 + .^1 - 2r2/r|) /2 = (19) 
= rs (^1 + s/l^^4S*J(N^^V^^^ /2, (20) 

where ([T9)) has been used. For a solar mass black hole, 
Ea.([T9l) implies that the number of charges needed to 
avoid the z — 1 singularity is just 

Nq = {2rsSl/lp)^^yV2^ ~ 2.91 x lO^^, (21) 



or 484 moles, which is a tiny amount on astrophysical 
terms. In general, Nq = Nq{M/Mq)'^/^ implies that in 
astrophysical scenarios ~ rs- 



Horizon location as A'^ changes 




Figure 1. The location of the external horizon is given by the 
intersection of the curve /i (solid green) with /2 (solid red 
curves labeled by Nq). The dotted curves correspond to GR, 
fi^ in blue and the various /l'^ in orange. Color online only. 

To study how the structure of these nonsingular black 
holes changes with size, it is useful to express 82 as 

= (22) 

Nq V aem 

According to (|10|) the horizons are located at the points 
where the curves 

h{z)^l + 5lG{z) (23) 

and 



/2 = z6lyj2a-/aemlNq (24) 

meet. In GR these curves are just fi^{z) — 1 — ^^/z and 
f2^{z) — ]2{z)l ^WZ. In Figdlwe see that the cut points 
for the external horizon are almost indistinguishable from 
those of GR for A'g > 30, while the inner horizon is ab- 
sent. Moreover, from the analysis of these curves, it is 
found that when > ^i, i.e., as Nq drops below the crit- 
ical value N^ = ^j2/aem ~ 16.55 (see Eq.jSSl)), the ex- 
ternal horizon disappears and the core becomes directly 
observable. From (fT9|) . the mass of these objects is 

M = Nl'^{2a,,nf'^mp/{^51) « Nl'^mp/^b, (25) 
and for the particular case we find 

3 /2 

= mp/{V25l) = ""fj « 1.23605mp. (26) 
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This number is also related with the total energy of 
the electric field in Born-Infeld nonlinear electrodynam- 
ics (see e.g. ,13,]). Here it manifests the regularizing 
role played by gravitation, which forces the electric field 
to remain in the z > 1 region. This suggests that 
the theory ([Ij may provide new insights on the prob- 
lem of sources in electrodynamics coupled to gravitation 
12| , which already arises in the Reissner- Nordstrom and 
Schwarzschild solutions of GR It is worth noting 

that if Nq is seen as an integer number, besides the core 
area, the mass of these nonsingular naked cores and black 
holes is quantized, as has been also recently claimed in 
[l5| . One thus expects a discrete spectrum of Hawking 
radiation, because physically allowed transitions should 
occur only between regular configurations. This illus- 
trates how Planck-scale physics may affect the perturba- 
tive predictions of the semiclassical approach. 



CONCLUSIONS AND PERSPECTIVES 

We have shown with an exactly solvable model that 
nonperturbative effects at the Planck scale can gener- 
ate important modifications on the innermost structure 
of black holes while keeping their macroscopic features 
essentially unchanged. The existence of constraints re- 
lating the scales of the theory (total mass and charge 
in our model) to provide nonsingular configurations, and 
their relation with the quantization of the geometrical 
structures involved, such as the core surface area, sug- 
gests that quantum black holes may have a much richer 
structure than their classical limits. The impact that our 
results could have for the experimental search of micro- 
scopic black holes [3l and naked cores makes it necessary 
a more complete analysis of the influence of the sources 
on the structure of these compact objects. This and other 
related studies are currently underway [12 1. 
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